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Abstract. We discuss how lensing by magnetic fields in galaxy clusters affects 
ultrahigh energy cosmic ray (UHECR) observations. As specific example, we use Virgo 
together with the cluster magnetic fields obtained earlier in a constrained simulation 
of structure formation including MHD processes. We find that, if M87 is the single 
source of UHECRs from Virgo, the emitted fiux is strongly anisotropic in the most 
interesting energy range, (50-100) EeV, and differs from the average value by a factor 
five or more for a significant fraction of observers. Since magnetic lensing is energy 
dependent, the external energy spectrum as seen by different observers varies strongly 
too. These anisotropics arc averaged out in the case that all active galactic nuclei 
in Virgo emit UHECRs. In both cases, the anisotropics of the emitted UHECR flux 
may introduce an important bias in the interpretation of UHECR data like, e.g., the 
determination of the source density Ug and the source energy spectrum of UHECRs. 
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1. Introduction 

The identification of tlie sources of ultrahigh energy cosmic rays (UHECRs) is one 
of the most longstanding problems of astrophysics. Only recently, evidence has been 
accumulating that we are at the dawn of "charged particle astronomy." A first piece 
of evidence have been the observation of the Greisen-Zatsepin-Kuzmin cutoff [1] in the 
UHECR spectrum with 5a significance by the HiRes experiment [2] which was recently 
confirmed by the AUGER experiment with 6a significance [3]. Secondly, anisotropies 
on medium scales have been found previously combining all available data of "old" 
CR experiments [1]. The data from the Pierre Auger Observatory (AUGER) presented 
in Ref. [5] confirmed these findings, showing also a surplus of clustering compared to 
an isotropic distribution in the broad range from 7 to 30 degrees. Such anisotropies 
were predicted as a consequence of the observed large-scale structure of matter [6] 
and favor therefore together with the presence of the Greisen-Zatsepin-Kuzmin cutoff 
the extragalactic origin of UHECRs. Comparing the expected angular auto-correlation 
function of different sources with the data, the authors of Ref. |7j suggested that the 
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sources of UHECRs are active galactic nuclei (AGN) or another strongly clustered sub- 
sample of galaxies. Finally, the Auger collaboration reported evidence for an anisotropy 
of the UHECR arrival directions, observing a correlation of the arrival directions 
of UHECRs with the positions of AGN or similarly clustered matter in the nearby 
Universe [8] . At present the significance of this correlation is only a the Sex confidence 
level and is not confirmed by the data from the HiRes experiment [9]. 

Most (auto-) correlation analyses assume identical sources, emitting isotropically 
UHECRs. Here we show that these assumptions are not fulfilled in the case that UHECR 
sources are located inside regions with relatively strong magnetic fields as in the core 
of galaxy clusters. As a concrete example we use the Virgo cluster and investigate how 
lensing effects of magnetic fields extending on cluster scales influence the emitted flux 
and energy spectrum both for isotropic and jet-like emission of single sources. After 
a brief discussion of cluster magnetic fields in Sec. 2, we examine in Sec. 3 first the 
anisotropies in the case of a single isotropic source, then of jet-like emission by a single 
source, of several isotropic sources and finally the resulting modifications of the observed 
CR energy spectrum. We conclude in Sec. 4. 

2. Magnetic fields in galaxy clusters 

Magnetic fields have been detected in galaxy clusters by radio observations, via the 
Faraday rotation signal of the magnetized cluster atmosphere towards polarized radio 
sources in or behind clusters [10] and from diffuse synchrotron emission of the cluster 
atmosphere (see Ref. [TTl [T2] for recent reviews). Although our understanding of their 
origin is still limited, at present, numerical simulations which follow the amplification 
of weak magnetic seed fields by structure formation lead to a consistent picture of the 
predicted magnetic field and its structure within galaxy clusters [I3l [lH [161 HZJi well 
in line with the observed Faraday rotation signal. It is worth mentioning that the 
predicted magnetic field structure within galaxy clusters is found to not depend on 
details within the structures of the magnetic seed fields and therefore can be assumed 
to be independent on the exact mechanism which creates the magnetic seed field at high 
redshift [HI [15]. 

To obtain the predicted magnetic field structure for the Virgo cluster we used 
the results from a constrained realization of the local universe (see [11] and references 
therein). In short, the initial conditions are similar to those used by Ref. [18] in their 
study of structure formation in the local universe. The initial density fluctuations were 
constructed from the IRAS 1.2-Jy galaxy survey by smoothing the observed galaxy 
density field on a scale of 7 Mpc, evolving it linearly back in time, and then using it as a 
Gaussian constraint [19] for an otherwise random realization of the ACDM cosmology. 
The volume constrained by the IRAS observations covers a sphere of radius ~ 115 Mpc 
centered on the Milky Way and therefore encloses the prominent local structure including 
the Virgo cluster. This region is sampled with high resolution dark matter particles 
and is embedded in a periodic box of around 343 Mpc length. The region outside the 
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constrained volume is filled with low resolution dark matter particles, allowing a good 
coverage of long range gravitational tidal forces. Many of the most prominent clusters 
observed locally can therefore be identified directly with halos in the simulation, and 
their positions and masses agree well with their simulated counterparts. For performing 
the MHD simulations, the original high resolution dark matter particles where split 
into gas and dark matter particles with masses of 0.69 x 10^ Mq and 4.4 x 10^ Mq, 
respectively. The gravitational force resolution (i.e. the comoving softening length) of 
the simulations was set to be 10 kpc, which is comparable to the inter-particle separation 
reached by the SPH particles in the dense centers of our simulated galaxy clusters. 

For this work we used the magnetic field configuration obtained from two 
realizations (MHDy and MHDz) starting from different initial seed fields (for more 
details see Ref. [Hj) having lower (MHDy) and higher (MHDz) initial values for the 
magnetic seed field. Both simulations lead to magnetic fields within the galaxy clusters 
which are statistical in agreement with the still rare Faraday rotation measures and are 
roughly bracketing the allowed range of magnetic fields in galaxy clusters. 

3. Anisotropies induced by cluster fields 

We consider first the case that only a relatively small subset of all AGN, i.e. of order 
1-10%, can accelerate to the highest energies, E > lO^'^eV. In particular, the number 
density of sources is then so small that in most clusters only one UHECR source is 
active. A concrete example for this situation is the acceleration of UHECRs in hot 
spots of radio galaxies [20]. For instance, there are only few radio galaxies inside the 
sphere with radius 70 Mpc used in the Auger analysis [8]. The two closest radio galaxies 
are Centaurus A and M87. The latter AGN sits in the center of the Virgo cluster. 

For the calculation of CR trajectories we use magnetic fields from the two MHD 
simulations of Ref. [H] . In these simulations, the magnetic field strength decreases fast 
for increasing distance from the cluster core, from fewxlO"^ G at the core over ~ 10~^ G 
at 0.5 Mpc to 10~^^ G at few Mpc. We use therefore a box of size 5 Mpc centered at the 
position of M87 and calculate the CR trajectories under the influence of the Lorentz 
force until the particles leave the box. Then we prolong the trajectory to obtain the 
final position on a sphere with radius d = 18 Mpc equal to the distance to the Earth, 
assuming that deflections outside the box can be neglected. This implies in particular 
that our results do not account for the effects of the Galactic magnetic field. We consider 
as two main cases first the isotropic emission of UHECRs and second jet-like emission. 

Before we discuss our results, we note that large deflections close to the source are 
not excluded by the Auger correlation signal. Since the extragalactic magnetic fields 
found in the simulations [H] are strongly concentrated in the core of the galaxy clusters 
on scales of size / ~ 0.5 Mpc, large deviations between the emission and exit direction 
from the cluster lead only to observed deflections § ^ l/d oi few degrees or less at the 
distance d = 18 Mpc to the Virgo cluster. More precisely, we find that the apparent 
size of a point source at the center of the Virgo cluster increases from 0.3° at 100 EeV 
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Figure 1. Skymaps of UHECR protons emitted uniformly from M87 witli energies 
E = 1000 (top, left) 100, 10, and 1 EeV (lower right panel). 

over 1.4° at 10 EeV to 3.6° at 1 EeV for an average observer, with sizeable deviations 
depending on the chosen direction. 

3.1. Isotropic emission of a single source 

We inject a large number of CRs with fixed energy E at the position of M87, choosing 
their initial propagation direction isotropically. Then we propagate each CR until it 
leaves a sphere of radius R = 2.5 Mpc around M87. We neglect all energy losses during 
the propagation, focusing only on the deflections due to the magnetic field. This is 
justified, because the average path-length in the box is even for iron at £^ = 10 EeV 
only 17 Mpc. Skymaps of the exit directions from this sphere are shown in Fig. [1] 
for protons with four different energies, E = 1000 EeV (upper left panel), 100 EeV, 
10 EeV, and 1 EeV (lower right panel), using the simulation with stronger magnetic 
fields (MHDz). Large anisotropics are clearly visible by eye in the most interesting 
energy range, (10-100) EeV. 

To quantify the probability to observe a fiux modified by magnetic lensing we use 
the following procedure: For any random exit directions around M87 we define small 
areas of size A ^ 7r5^ and calculate the number iVj of CRs per area and the average 
value (A^). Then we count the fraction f{e) of areas where the relative number Ni/ (N) 
of CRs is below or above a certain threshold value e. The function f{e) depends on the 
chosen resolution 6, i.e. on the size A = 716"^ of the areas used for averaging, /5(e). 
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Figure 2. The fraction of regions with flux above or below a certain threshold value 
for several different viewing angles 5 together with the extrapolation curves to (5 ^ 0; 
Red, solid line is for regions with flux twice smaller, blue dashed line for twice larger 
and magenta dotted line for five times smaller flux than average. 

A CR experiment collects during the course of a year CRs from Virgo within an 
opening angle d ~ 2AU/18Mpc ^ 10~^^ on a sphere centered in Virgo and radius 
18Mpc. We calculate therefore f&{s) for various values of 5 and extrapolate then to 
5^0. This requires a large enough total number of cosmic rays, so that all Ni 
are sufficiently high. For fixed 5, e and energy both the value of f5{s) and its 
uncertainty depend on the number of MC simulations. In order to extrapolate 5^0, 
we need to know not only /5(e), but also its fiuctuations due to the finite number of 
MC simulations. In order to reduce the computing time, we approximate the error as 
df = dN/Ntot = f I \J {Ni). By performing many MC simulations for several energies we 
checked that this approximation for df is a factor 2 or less within the exact one. 

An example for the extrapolation procedure 5 ^ is shown in Fig. [2], where the 
probability to observe a fiux two times smaller (solid line), two times larger (dashed line) 
and five times smaller (dotted line) than the average fiux is shown for model MHDy 
and E = 50 EeV for five different values of the areas A ^ 7r5^. The errorbars show the 
error due to the finite number of MC simulations as discussed above. Fitting a second 
order polynomial a + b6 + cS^ to these values, we obtain finally the asymptotic value of 
fb{e) = a for 5 — s> from a fit. 

The results for f{e) obtained in this way using the two different magnetic fields are 
shown in Fig. [3l Although the strength of the magnetic field at the cluster center and the 
deflection of individual CR trajectories differs considerably between the two simulations, 
the anisotropy of the emitted CR flux is very similar, in particular at low energies. The 
reason for this result becomes clear if one follows individual CR trajectories: The main 
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Figure 3. The fraction of sky /(e) exceeding the indicated level of anisotropy, for 
different energies and for magnetic fields from MHDy (left) and MHDz (right panel). 

part of the deflections is caused by the extended field outside the cluster core, that has 
roughly the same strength in both simulations. 

We discuss now the results obtained for different energies. For E = 1000 EeV 
protons propagate in the ballistic regime and deflections are small. As a result, the 
sykmap of the emitted CR flux shown in the upper left panel of Fig. [1] is essentially 
isotropic, while the fraction of sky with a flux modulated by a factor two is only on the 
1% level. The situation is already qualitatively different at the for us most interesting 
energy, E = 100 EeV. Deflections become sizable and regions with clearly reduced or 
enhanced flux are visible on the skymap, cf. the second panel of Fig. [H The CR flux is 
lower than average by a factor two in 30% and by a factor 5 in 5% of the sky. 

At E = 10 EeV protons start to propagate in the quasi-diffusive regime, and the 
anisotropies are weaker. However, broad over- and underdense regions are still visible in 
Fig. [TJ By contrast, at E = 1 EeV protons propagate in the truly diffusive regime and 
most anisotropies caused by the magnetic field component are washed out: Anisotropies 
on a level of a factor two occur still in 4% of the sky, while stronger anisotropies have 
such a small probability that we could not determine them reliably (/ < 10~^). 

The expected anisotropies for the case of iron (charge Q = 26) or other nuclei can 
be read off from Fig. [3] by the appropriate rescaling of their energy, E — E/Q. Thus 
strong anisotropies, e.g. larger than a factor five, are already unlikely for iron nuclei 
with E = 100 EeV and cluster magnetic fields of the strength as typically predicted 
by the simulations of Ref. [H]. On the other hand, the CR flux is also for iron nuclei 
modulated by a factor of two in a large fraction of the sky. 

We conclude that if M87 is the single source emitting UHECR in the Virgo cluster, 
then the emission as seen from outside of the cluster is rather anisotropic in the most 
interesting energy region both for iron and proton primaries. For a significant fraction of 
observers the flux differs from the average value by a factor flve (protons) or two (iron), 
respectively. Hence, in models assuming that only a small fraction of all AGNs can 
accelerate UHECR it would be natural to observe no UHECRs from a certain fraction 
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Figure 4. Skymap of UHECRs emitted in two beams (jet-like emission) from M87 
with energy 60EeV, red for protons and green for iron. The original directions are 
blue, the left panel uses the magnetic field from the MHDy and the right from the 
MHDz simulation. 

of galaxy clusters, as long as the total number of UHECR observed is as small as at 
present. In particular, it would be no contradiction to the AGN source hypothesis that 
no UHECRs from the Virgo cluster have been observed with the present statistics. 

3.2. Jet-like emission of a single source 

We consider now the case that the acceleration of cosmic rays results in the beamed, 
jet-like emission of UHECRs [21]. While in the previous case of isotropic emission the 
magnetic field created anisotropies, its effects is now opposite, washing out the strong 
initial anisotropy. We compare again the results obtained for two different magnetic 
field models from Ref. [H], in the left and the right panels of Fig. HI Each panel shows 
two jets with initial opening angle 5° in blue as well as the deflected exit directions 
in red for protons and in green for iron, both with E = 60EeV. Both the direction 
and the magnitude of the deflection varies between the two simulations and different 
sky regions. However, both simulations agree that for the energy region relevant for 
UHECR astronomy, E > Eqzk ~ 50 EeV, the solid angle Qe« of significant UHECR 
emission is increased by not more than a factor two for protons, while the beaming of 
iron nuclei is destroyed nearly completely, cf. Fig. HI Hence beaming will decrease the 
effective number density of iron sources only marginally, while this effect is prominent 
for proton sources. 

Since the clustering properties of the UHECR arrival directions [TJ [25] and the 
correlation analysis [8] of the Auger collaboration, taken at face value, favour that the 
fraction of all AGN visible to us is large, the possibility that only a small subset of all 
AGN like radio galaxies emits narrow proton beams is strongly disfavoured. 
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3.3. Many sources per cluster 

We assume now that all AGN emit UHECR and chose from the 12th VCV catalogue [22] 
all AGN closer than 2.5 Mpc to M87 as sources. Within this distance there are 
13 additional AGN and many of them are outside the core of the Virgo cluster in 
regions with low magnetic fields. Moreover, the distance between the AGNs is much 
larger than the correlation length of the field and therefore the over- and underdense 
regions produced by different AGNs are un-correlated. Consequently, the anisotropies 
are drastically reduced compared to the single source case. 

Assuming that all AGN have the same luminosity, we obtain a very uniform 
distribution of outgoing UHECR in the case of protons primaries: In less than 1% 
of all the sky the flux deviates more than a factor two from the average flux of UHECR 
protons with E = 100 EeV. Using iron nuclei with the same energy increases slightly 
the anisotropies: The iron flux is changed by a factor 2 in around 4% all of areas. We 
stress however that anisotropies induced by magnetic lensing for single sources lead to an 
important bias in studies of the density of UHECR sources and their energy spectrum, 
although the flux emitted by the total cluster is rather isotropic. 

3.4- Modified external CR energy spectra 

Since the (de-) magniflcation of the CR flux by magnetic flelds is energy dependent, the 
anisotropies observed in the previous subsections lead to a modulation of the initial CR 
spectrum for an outside observer. In the Fig. [5] we show the energy spectrum dN/ dE 
as function of energy for 10^ CRs and three different regions of size 5 = 8° in the 
top panels, while skymaps of 5000 UHECR protons emitted with an energy spectrum 
dN/dE oc i?"^'^ and threshold energy E > Eq = 5 EeV are shown in the bottom. We 
assume flrst that only one source, M87, emits UHECRs. For both models MHDy (left) 
and MHDz (right) we show the original spectrum in red, in a region with average flux 
(blue crosses), and in two regions with increased (green boxes) and reduced (magenta 
circles) flux, respectively. In all cases, the spectrum is larger or smaller than average 
for several nearby bins. This mean that lensing effects are not washed out by the 
finite experimental energy resolution. Also the statistical errors shown in Fig. [5] are 
small enough not to affect the conclusions above. Thus, magnetic lensing does not 
only change the value of the UHECR flux in a given direction, but also changes its 
energy dependence: Depending on the threshold energy and the direction considered, 
the exponent of the estimated power-law dN/dE oc i?"" may differ strongly. Note that 
the shown energy spectrum will be additionally modifled both by propagation effects 
and by the Galactic magnetic fleld. The latter will not only deflect CRs but will lead 
additionally to a similar modulation of the CR energy spectra as the cluster fleld. 

Next we investigate how fast this modiflcation of the CR spectrum disappears, if the 
the number of sources increases. In Fig. [6l we show again the energy spectrum for two 
regions with reduced (left) and increased (right) flux, varying now the number of sources. 
While in the case of a spot with small anisotropy the modifled energy spectrum resembles 
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Figure 5. The energy spectrum dN/dE as function of energy (top) and skymaps 
of UHECR protons emitted uniformly from M87 (bottom) for models MHDy (left) 
and MHDz (right). The energy spectrum is shown for 5 = 8° regions and 10^ events 
over the sky for the following cases: in absence of magnetic field (red), spectrum in 
overdense region (magenta), underdense region (green) and ordinary region (blue). 

already for two sources the original one, the modification of the energy spectrum in case 
of a large anisotropy is non-negligible even if all 14 AGN contribute equally to the 
flux. Finally, we remark that the case with 14 sources with equal luminosity is rather 
unrealistic, since the large difference in (bolometric) luminosity between M87 and the 
other AGN in the Virgo cluster should be reflected also in their UHECR luminosity. 

4. Conclusions 

The amplification of the UHECR flux by magnetic lensing was discussed previously in 
the context of the Galactic regular and turbulent field [23l[2l]. For instance, Ref. [23] 
calculated exposure maps for three different models of the Galactic regular magnetic field 
and found also significant anisotropies at = 40 EeV. The strength of the magnetic 
field close to the core of galaxy clusters is smaller by a factor 30-100 than in the Milky 
Way. The weakness of cluster fields is, however, over-compensated by their much 
larger extension compared to the Galactic disc. While the deflection of UHECRs in 
cluster fields leads even for the nearest cluster only to maximal angular differences 
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Figure 6. The energy spectrum dN/dE as function of energy for model MHDz and 
a varying number of sources; left for underdense and right for an ovcrdcnsc region. 

between the source and the arrival direction of order 1.5° and is thus neghgible, the 
resulting anisotropics of the emitted CR flux may introduce an important bias in the 
interpretation of CR data. 

We have found that assuming magnetic fields as in Ref. [Hj the flux from M87 
differs from the average value by a factor five for a significant fraction of observers in the 
most interesting energy region, (50-100) EeV. This may be the explanation for the non- 
observation of UHECR from the Virgo cluster, if M87 is the dominant UHECR source 
in the Virgo cluster. More generally, anisotropics induced by magnetic fields result in a 
reduction of the effective number of sources, leading to a bias in auto-correlation studies. 
Vice versa, if the source type of UHECRs and thus their number density is known, a 
comparison with the estimated value from auto-correlation studies informs us about 
the importance of anisotropics and/or beaming effects. Note that magnetic lensing will 
affect also sources behind the Virgo cluster. 

The second important consequence of the anisotropics induced by the cluster 
magnetic fields is the resulting modulation of the energy spectrum of UHECRs leaving 
the galaxy cluster compared to their generation spectrum. A similar energy modulation 
is be introduced by the Galactic magnetic field and both effects will complicate the 
reconstruction of generation spectrum, if UHECRs sources are discovered. 

Finally we want to stress that we obtained all our results for two specific realizations 
of magnetic fields obtained within the simulations of Ref. [H]. Although we expect that 
the main features of our results are robust, quantitative differences will appear using 
different simulations of cluster magnetic fields: For instance, if the volume fraction filled 
by relatively large magnetic fields around clusters is as large as in Ref. [26], we expect 
that the same phenomena happen as found in this work, but extending over a larger 
energy range. 
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